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1. Introduction 

The concept of pseudospin symmetry [1, 2] is based on the empirical ob- 
servation of quasi-degenerate pairs of certain normal-parity shell-model or- 
bitals with non-relativistic quantum numbers 

nr,t,j = £ + ^) and (n r - 1, I + 2, j = I + |) . (1) 

Here n r , £, and j are the single- nucleon radial, orbital, and total angular 
momentum quantum numbers, respectively. The doublet structure, is ex- 
pressed in terms of a "pseudo" orbital angular momentum I = £ + 1 and 
"pseudo" spin, s = 1/2, coupled to j = £±s. For example, (2p 3 / 2 , 1/5/2) will 
have £ = 2, etc. This pseudospin symmetry plays a central role in nuclei [3] 
and only recently has it been shown to originate from a relativistic symme- 
try of the Dirac Hamiltonian [4, 5]. The members of the pseudospin doublet 
exhibit the following features, (i) They have different angular momentum 
quantum numbers (£, j) <-» (£ + 2, j + 1) for states with aligned/unaligned 
spin, (ii) They have different radial wave functions and, in particular, dif- 
ferent number of nodes n r <-> n r — 1. (iii) They involve only normal-parity 
shell- model orbitals. The "intruder" levels with aligned spin and no nodes, 
e.g. Ogg/2, 0/in/ 2 , 0« 13/ / 2 do not form quasi-degenerate doublets. In the 
present contribution we show that a natural explanation for all these fea- 
tures can be obtained by combining the relativistic attributes of pseudospin 
symmetry with known properties of Dirac bound states [6] . 
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2. Dirac Wave Functions in Central Fields 

The eigenstates of a spherically symmetric Dirac Hamiltonian have the form 
* K ,m = [Y £ x] ( £ , i^f [Ye x}$ ) where G K {r) and F K (r) are the radial 
wave functions of the upper and lower components respectively. The labels 
« = T(j + 1/2) and £' = I ± 1 for j = I ± 1/2. As bound state solutions, 
both G(r) and .F(r) vanish at r = and r = oo and the binding energy is 
positive (M — E 1 ) > 0. Further properties of these wave functions can be 
inferred from the following equation for the product GF 

(GF)' = A{r) F 2 - B(r) G 2 , (2) 

where A(r) = [E + M + V s {r) - V v {r)], B(r) = [E - M - V s (r) - V v {r)]. 
For relativistic mean fields relevant to nuclei, the scalar potential, Vs(r), is 
attractive and the vector potential, Vv(r), is repulsive. For typical values 
A(r) > 0, while B(r) changes monotonically from -B(O) > to B(oo) = 
{E — M) < 0. Since GF vanishes at both end points, we find by Eq. (2) 
that it is an increasing negative function at large r, while at small r, GF is 
a decreasing negative function for k < and an increasing positive function 
for k > 0. Futhermore, since B(r) > at the nodes of F and G (nodes can 
occur only in the region where the kinetic energy is positive), it follows that 
GF is a decreasing function at the nodes of F, and an increasing function 
at the nodes of G. Exploiting all these properties, we observe that for k > 0, 
GF is positive at small r and negative at large r, and hence has an odd 
number of zeroes. The first and last zeroes of GF correspond to nodes of F 
and since the nodes of F and G alternate [6, 7], then the number of nodes 
of F exceed by one the number of nodes of G. On the other hand, for k < 0, 
GF has the same (negative) sign near both end points, and hence has an 
even number of zeroes. By similar arguments we find that in this case the 
first and last zeroes of GF are nodes of G and F respectively, and that G 
and F have the same number of nodes. These properties of the radial wave 
functions are illustrated in Figs, la-d for the 2p 3/ / 2 (k = —2), I/5/2 (k = 3) 
states. Altogether we have, 

k < : n F = n G , . 

k > : np = no + 1 

where np and ng denote the number of internal nodes of F and G. 

When the radial parts of the wave functions G and F have no nodes 
(n.F = no = 0), the corresponding bound states can appear only in the 
j = i + 1/2 state (k < 0), but not in the j = i - 1/2 state (n > 0) [6, 8]. 
The reason being that in this case GF has the same sign for all r and since 
always GF < at large r, then it is also negative at small r and therefore 
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Figure 1. (a) The radial upper component (G) and lower component (F) in (Fermi) -1 / 2 
and (b) the corresponding product GF in (Fermi) -1 of the 2p 3 / 2 [k = —2) state, (c) and 
(d) The same for the I/5/2 {k = 3) state, (e) and (f) The same for the 0g a / 2 ( K — — 5) 
state. All states are eigenstates of a Dirac Hamiltonian with scalar (S) and vector (V) 

potentials: Vs,vi r ) = a s,v [l + exp(^5)] 1 with parameters: as = —358, av = 292 
MeV, R = 7 fm, a = 0.6 fm, tuned to the neutron spectra of 208 Pb. Taken from Ref. [6]. 
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must have n < 0. An example is shown in Figs, le-f for the case of the 0g$/2 
(n = —5) intruder state. 

Finally, we note that for bound states to exist, there must be a region 
where B(r) > 0. This is due to the fact that, as mentioned, GF is an 
increasing negative function at large r, and to enable it to vanish at r = 0, 
its derivative (GF)' must change sign from positive to negative in some 
region. A glance at Eq. (2) shows that since A(r) > 0, a necessary (but not 
sufficient) condition for (GF)' to become negative is that 

B(r) = [E-M - V s (r) - V v (r) ] > for some r . (4) 

The above condition means that in order that bound states exist, there 
has to be a region where the depth of the average attractive single-nucleon 
potential, Vs(r) + Vy (r), is larger than the binding energy. 



3. Relativistic Pseudospin Symmetry 

All of the above results are relevant for understanding properties of states 
in the relativistic pseudospin scheme [4, 5]. When Vs = —Vy, the Dirac 

Hamiltonian is a scalar under an SU(2) algebra, [H , S M ] = 0. If in addi- 
tion the potentials are spherically symmetric, the Dirac Hamiltonian has 

an additional invariant SU(2) algebra, [ H , ] =0. The relativistic pseu- 
dospin generators, S^, and relativistic pseudo-orbital angular momentum 
operators, L M , are given by 

Mu)<MH)' 

Here = U p U p and = U p U p , where = and 1^= rx p are 
the usual spin and orbital angular momentum operators respectively, the 
Pauli matrices and U p = — is the momentum-helicity unitary operator 

introduced in [9]. The sets js^, s M , s M | and i^, j form two triads 

of SU(2) algebras. The and operators act on the four-components 

Dirac wave functions. The s M and operators form the non-relativistic 
pseudospin and pseudo-orbital angular momentum algebras respectively, 
and act on the upper components of the Dirac wave functions. The s M and 
act on the "small" lower components of the Dirac wave functions. The 
pseudospin s = 1/2 and pseudo-orbital angular momentum I are seen from 
Eq. (5) to be the ordinary spin and ordinary orbital angular momentum 
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of the lower component of the Dirac wave functions. The two states in the 
pseudospin doublet share a common pseudo-orbital angular momentum £, 
which is coupled to a pseudospin s = 1/2, and thus have the form 

*«<„,„ = (^2W K _ lX ]«,^ Kx ]«) j = i-\ 

*V>o, m = {^^[Y i+ ,X&\'^^[Y i x]t ) )f = i+\- (6) 

The underlying Dirac structure ensures that the wave function of the up- 
per component of the Dirac eigenfunction has a spherical harmonic of rank 
£ = £ — 1 for aligned spin: j = £ — 1/2 = £ + 1/2, and a spherical harmonic 
of rank £ + 2 = £ + 1 for unaligned spin: / = £ + 1/2 = {£ + 2) - 1/2. This 
explains the particular angular momenta defining the pseudospin doublet 
of Eq. (1). In the pseudospin symmetry limit the two states in Eq. (6) 
form a degenerate doublet (5 = 1/2), and are connected by the pseudospin 

generators 5 M of Eq. (5). The corresponding upper components are a dou- 
blet with respect to the set (the non-relativistic pseudospin algebra). 
Since the latter, by definition, intertwine space and spin, they can connect 
states for which the upper components have different radial wave functions, 
G K <o(r) ^ Gv>o(r). On the other hand, the corresponding lower compo- 
nents are a doublet with respect to the ordinary spin s^, and hence, in the 
pseudospin limit, their radial wave functions are equal up to a phase, 

F K<0 (r) = F K/>0 (r) . (7) 

In particular, F K< o(r) an d -?K'>o( r ) have the same number of nodes, which 
we denote by n r . If we now use the result of Eq. (3), we find for n r ^ that 
G K <o( r ) in Eq. (6) has also n r radial nodes, while G K />o(r) has n r — 1 nodes. 
This explains the structure of nodes in the pseudospin doublet of Eq. (1). 
The simple relation in Eq. (7) between the radial wave functions of the 
lower components of the two states in the doublet, dictates this particular 
relation between the radial nodes of the corresponding upper components. 
This result cannot be obtained if one considers just the non-relativistic 
pseudospin algebra, and is a direct outcome of the behavior of nodes 
of Dirac bound states and the identification of pseudospin as a relativistic 
symmetry of the Dirac Hamiltonian. 

As we have shown, bound Dirac states, for which both the upper and 
lower components have no nodes (n r = 0) can occur only for k < and not 
for k' > 0. From Eq. (6) we find that such states have pseudo-orbital angular 
momentum £ and total angular momentum j = £—1/2. As mentioned, these 
intruder states are ignored in the non-relativistic pseudospin scheme, and 
it is only the relativistic interpretation of pseudospin symmetry, combined 
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with known properties of Dirac bound states, which enable a classification 
for these states, as well as provide a natural explanation why these states 
do not have a pseudospin partner which is an eigenstate of the Hamiltonian. 

The exact pseudospin limit requires that Vs(r) = —Vy(r), which implies 
that B(r) = E — M. It is clear that under such circumstances the condition 
of Eq. (4) cannot be fulfilled for bound states with positive binding energy 
M — E > 0. This explains why in the exact pseudospin limit, there are 
no bound Dirac states and, therefore, by necessity the pseudospin symme- 
try must be broken in nuclei. Nevertheless, a variety of realistic mean field 
calculations show that the required breaking of pseudospin symmetry in 
nuclei is small [10, 11, 12, 13]. Quasi-degenerate doublets of normal-parity 
states, and abnormal-parity levels without a partner eigenstate persist in 
the spectra, and the relation of Eq. (7) is obeyed to a good approxima- 
tion, especially for doublets near the Fermi surface. As discussed, these 
features are sufficient to ensure the observed structure of nodes occurring 
in pseudospin doublets and the special status of intruder levels in nuclei. 

4. Acknowledgements 

It is a pleasure to dedicate this article to P. Ring on the occasion of his 60th 
birthday. This research was supported in part by the U.S. -Israel Binational 
Science Foundation and in part by the United States Department of Energy 
under contract W-7405-ENG-36. 

References 

1. A. Arima, M. Harvey and K. Shimizu, Phys. Lett. B 30 (1969) 517. 

2. K.T. Hecht and A. Adler, Nucl. Phys. A 137 (1969) 129. 

3. J.N. Ginocchio and A. Leviatan, contribution to these proceedings and references 
therein. 

4. J.N. Ginocchio, Phys. Rev. Lett. 78 (1997) 436. 

5. J.N. Ginocchio and A. Leviatan, Phys. Lett. B 425 (1998) 1. 

6. A. Leviatan and J.N. Ginocchio, submitted to Phys. Lett. B, nucl-th/0108016. 

7. M.E. Rose and R.R. Newton, Phys. Rev. 82 (1951) 470. 

8. M. Hirooka, G. Konisi and T. Satto, Prog. Theor. Phys. 41 (1969) 161. 

9. A. L. Blokhin, C. Bahri, and J.P. Draayer, Phys. Rev. Lett. 74 (1995) 4149. 

10. J.N. Ginocchio and D. G. Madland, Phys. Rev. C 57 (1998) 1167. 

11. GA. Lalazissis, Y.K. Gambhir, J.P. Maharana, C.S. Warke and P. Ring, Phys. Rev. 
C 58 (1998) R45. 

12. J. Meng, K. Sugawara-Tanabe, S. Yamaji, P. Ring and A. Arima, Phys. Rev. C 58 
(1998) R628; J. Meng, K. Sugawara-Tanabe, S. Yamaji and A. Arima, Phys. Rev. 
C 59 (1999) 154. 

13. J.N. Ginocchio and A. Leviatan, Phys. Rev. Lett. 87 (2001) 072502. 



